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ABSTRACT 

When an isolated magnetar with magnetic dipole field B ~ 10 15 G moves at high velocity (v > 10 7 
cms -1 ) through the ISM, its transition into propeller effect driven spin-down may occur in less than 
10 6 years. We propose that the nearby neutron star RX J1856. 5-3754 is such a magnetar, and has spun 
down by the propeller effect to a period greater than 10 4 sec within ~ 5 x 10 5 years. This magnetar 
scenario is consistent with observed thermal X-ray emission properties and the absence of detectable spin- 
modulations of them. Detection of other rapidly moving long period (> 100 sec) magnetars with known 
ages would strongly constrain the very great variety of predicted propeller-effect torque magnitudes. 

Subject headings: stars: neutron - X-rays: stars - stars: RX J1856. 5-3754 



1. INTRODUCTION 

Isolated magnetars (spinning-down neutron stars with 
dipole fields 10 14 -10 15 G) differ from canonical (10 n -10 13 
G) pulsars in several ways: 

(1) they spin-down very much more rapidly; 

(2) when relatively young (age 10 6 years), their spin- 
down torque mechanism may change over from the rela- 
tivistic particle plus magnetic field wind of canonical pul- 
sars to a "propeller" interaction (Davidson & Ostriker 
1973; Illarionov & Sunyaev 1975) on a surrounding inter- 
stellar medium if they are rapidly moving through it; 

(3) their strongly magnetized surface can have very dif- 
ferent emissivities for extraordinary and ordinary mode 
thermal X-rays. 

We consider below the extraordinarily varied spin-down 
descriptions for magnetars when various proposed pro- 
peller interaction models are used (Illarionov & Sunyaev 
1975; Davies et al. 1979). Especially significant may 
be a possible almost complete quenching of the spin of 
some magnetars within 10 years, not possible for iso- 
lated canonical pulsars or for magnetars embedded in their 
own accretion disks. We apply such considerations to un- 
derstanding the observed X-ray emission from the rapidly 
moving isolated 10 6 year old neutron star RXJ1856. 5-3754 
(hereafter RXJ1856). It has featureless nearly blackbody 
thermal X-ray emission (with a paradoxically small black- 
body radius) which shows no indication of spin-period 
modulation (< 1 %) in the range 10 ms < P < 10 4 sec 
(Ransom et al. 2002; Drake et al. 2002; Burwitz et al. 
2003). 

Four possible explanations of the lack of detected spin- 
modulation of the surface X-ray emission of an isolated 
10 6 year old NS are 

(1) Our line of sight to this pulsar is almost exactly along 
the pulsar spin-axis (Braje & Romani 2002); 

(2) the magnetic field of this neutron star, unlike that 
of all known pulsars, is almost exactly axially symmetric 
with respect to the neutron stars spin-axis (Braje & Ro- 
mani 2002); 

(3) this particular neutron star was born less than 10 6 
years ago as a weak field millisecond pulsar (P < 10 ms) 

1 Columbia Astrophysics Laboratory, New York, NY 10027 

2 Columbia University, Physics Department, New York, NY 10027 



(4) this neutron star's spin has been almost quenched 
to P > 10 sec. Of these (4) may be the most plausible. 
Whatever the details of the mechanism for achieving it are, 
the needed effective spin-down torque at very long periods 
seems to require both a transition out of the canonical pul- 
sar spin-down mechanism, where energy and angular mo- 
mentum are lost from an almost corotating magnetosphere 
at its very distant ( ^ 10 14 cm) light cylinder radius, and 
a huge surface dipole magnetic field. These two require- 
ments encourage exploring the possibility that RX J1856 is 
a magnetar, rapidly moving through the ISM or a molec- 
ular cloud, is now spinning-down by propeller ejection of 
that incident matter. 

The propeller effect may be an efficient spin-down mech- 
anism when the spin-period becomes larger than seconds, 
but proposed spin-down torques from it are highly model- 
dependent. We find that RX J1856 should spin-down to a 
period > 10 6 sec for some proposed propeller-effect mod- 
els. The magnetar scenario is also consistent with observed 
spectral properties for both a single temperature and two- 
component blackbody model. 

2. SPIN-DOWN EVOLUTION 

When an isolated neutron star (NS) is spun-down by 
the torque of its relativistic wind emission, 

m~-!t^L = -KSi n , (i) 

where / is the NS moment of inertia, fi is the NS magnetic 
dipole moment and Q is the NS spin-velocity. The "brak- 
ing index" n is 3 if /x and fi ■ Q,/Cl are constant. When 
the NS is born with a spin-period much shorter than a 
presently observed one, the age of the star is the "spin- 
down" age T sc i = 7— T7^, where P and P are the present 
spin period and its time derivative. The canonical spin- 
down age, r cs d = P/2P, differs from the true age for n^3. 
If rc < 1 r s d clearly must have a very different meaning. 
When 7i > 1, the canonical spin-down age indicates the 
time since the NS entered the spin- down phase described 
by equation (1): 

2 

T~sd = TTcsd, (2) 
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as long as the original spin-period <C P but when n < 1, 
the canonical spin-down age indicates the approximate 
time remaining before the NS spin will be completely 
quenched (rf): 

Tf = — T cs(i . (3) 

To quench a NS's spin, that spin-down must then be in a 
phase described by n < 1 in equation (1). 

When the NS light cylinder radius (c/O) becomes larger 
than the radius of its corotating magnetosphere (R m ), the 
spin-down mechanism is changed to one in which the spin- 
down torque is from direct spin-up of the surrounding 
medium ("propeller effect"). In this case, the braking in- 
dex is no longer ~ 3. In proposed models, it is generally 
in the range — 1 < n < 2. 

Propeller effect spin-down torque depends upon three 
parameters, the NS magnetic dipole moment (p), the den- 
sity of incoming gas at the magnetosphere boundary with 
the surrounding medium (p m ) and velocity of the incom- 
ing gas at the magnetosphere boundary in the NS rest 
frame (v m ). Assuming the NS spin-down torque has fixed 
power-law dependences on these parameters, dimensional 
analysis gives 

m = — k ^ +n)lz pt- n)l% vt^ n)lz si n (4) 

with k and n dimensionless constants. In the following 
discussion, we set NS mass and radius to 1.4M and 10 
km respectively. 

2.1. Models for propeller caused spin-down 

Two parameters can classify the many proposed pro- 
peller spin-down models. In these models, spin-down 
torque is 

ISl = -MR m v m M\ (5) 

where M = Pm(^Rm) v m is the incoming (and largely 
ejected) mass rate. M. is the Mach number defined as the 
ratio of incoming medium velocity to NS spin-velocity at 
the magnetosphere boundary: M. = R m Q/v m . Proposed 
propeller models have 7 between —1 and 2. M > 1 and 
M. < 1 are the supersonic and subsonic regimes, which 
have been treated differently in proposed models. 

The ram pressure balance at the (non-spherical) mag- 
netosphere boundary from either v m or £!i? m is assumed 
to be 
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where B m is local magnetic field strength at the magneto- 
sphere boundary. Since B m — pjR? m for a dipole magnetic 
field with dipole moment p, the braking index n is then 

3(2 7 -<S) 
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The spin-down torque is 
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A remarkably large number of different propeller spin- 
down models have been proposed and applied. The brak- 
ing indices suggested by these propeller models are in the 
range of — 1 < n < 2. Some of the models were intended 
only for NSs with spin-down from surrounding accretion 



disks. Table 1 categorizes most of the proposed propeller 
effect models by their 7 and 5. 

We note that at the transition into propeller spin-down 
of an isolated pulsar, f2i? m ~cso that M <~ 10 3 in equa- 
tions (5), (6) and (8). In the models of Table 1 (by vari- 
ous authors, and, sometimes, the same authors at different 
time) the j-5 exponent in equation (8) varies from —1 to 2 
corresponding to differences of 10 9 in the magnitude of the 
propeller driven spin-down torques which take over from 
that of equation (1). 

In most models 7 7^ S so that in the transition from 
pulsar to propeller phase there is a huge discontinuity in 
spin-down torque. We favor a plausible propeller spin- 
down model with 7 = 5 which gives no such discontinuity. 

A particularly simple model completely ignores all mag- 
netic fields within both the incoming and ejected ionized 
beams and also possible collisionless collection stream- 
stream interactions. Because both flows are so very dilute 
{n m ^ 10 cm -3 ) single particle interaction mean free paths 
3> i£ m (~ 10 12 cm). Incidence and reflection of the ISM on 
the stellar magnetosphere could then be described as a sum 
of single particle trajectories. A "rough" interface model 
(RIM) is one in which the interface between the corotat- 
ing magnetosphere of the NS and the ISM flow onto it is 
approximated as a sawtooth-like structure in the corotat- 
ing frame, with angles between the sawtooth surfaces and 
the surface of constant radius from the star assumed to 
be of order a radian. In "smooth" interface models that 
interface is assumed to be much more spherical (typically 
with angular deviations of order the inverse Mach number 
of the very supersonic propeller motion) . In a " rough" in- 
terface model the deflected incident flow pushes strongly 
with a comparable force in the tangential and radial di- 
rections. This gives 7 = 5 = 1 and n = 3/7. (In a smooth 
interface limit, 7 = — 1,5 = gives the n = — 1 of Illari- 
onov & Sunyaev (1975)). In our discussions below we give 
particular emphasis to it. 

2.2. Transitions between different spin-down stages of 
magnetars 

We discuss three transitions in the spin evolution of 
magnetars as ft decreases. When ft = c/R m there is a 
transition from canonical pulsar spin-down to propeller 
spin-down. When ft = v m /R m , there is the transition 
from supersonic (M > 1) to subsonic propeller spin-down 
(A4 < 1). When £1 = f2^(i? TO ) (the Keplerian angular 
velocity at the magnetosphere boundary), there is a tran- 
sition from propeller spin-down to strong accretion. Table 
2 summarizes spin-periods and transition times for differ- 
ent propeller spin-down models. 

3. THE EXCEPTIONAL ISOLATED NEUTRON STAR 

RX J1856 

The isolated neutron star RX J 1856 was serendipitously 
discovered by ROSAT as an X-ray source with thermal 
emission (Walter et al. 1996). Detection of proper motion 
provided a distance ~ 120 pc and velocity <~ 200 kms -1 
and supports an association of RX J1856 with the Upper 
Scorpius birthplace indicating a NS age <~ 5 x 10 5 yrs (Wal- 
ter & Lattimer 2002; Kaplan ct al. 2002). Deep Chandra 
and XMM-Newton observations with total exposure time 
~ 500 ks were performed to search for spectral feature and 
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pulsation. However, neither was found in the X-ray data 
which had high resolution and sensitivity (Ransom et al. 
2002; Drake et al. 2002). 

This lack of observed pulsation is perplexing since mod- 
ulation of an anisotropic temperature distribution sug- 
gested by the two-component blackbody model fit dis- 
cussed below seems expected, and pulsation modulation 
of X-ray emission has been detected in other isolated NSs 
with thermal emission (Haberl et al. 1997; Zavlin et al. 
2000). We propose that RX J1856 has been spun-down to 
a spin period longer than 10 4 sec by the propeller effect. 
The fraction of magnetars among isolated NSs is ~ 10~ 4 
(Kouveliotou et al. 1998), and is much larger among NSs 
from which thermal emission has been detected. We find 
that the only plausible way to achieve the needed condi- 
tion for early transition into the propeller phase and rapid 
enough spin-down thereafter is for this star to be a mag- 
net ar. 

We first assume B > 10 14 G, sufficiently large that tran- 
sition into the propeller phase occurs in much less than 
5 x 10 5 years. Then, we searched parameter space in the 
j-5 plane for the region in which RXJ1856 could, there- 
after, have been spun-down to a period longer than 10 4 
sec. Input parameters are magnetic dipole moment /i, in- 
coming gas density n m , and incoming gas velocity v m . We 
fixed v m from the observed stellar proper motion at 200 
km _1 . 

ISM particles (mainly hydrogen) at the magnctosphere 
radius will be charged and be reflected by the stellar mag- 
netosphere. [Thermal photons from the NS surface ionize 
all hydrogen in the vicinity of the star. The ionization 
time of hydrogen at R m <~ 10 12 cm (<~ 10 2 s) is much 
shorter than R m /v m (~ 10 5 s). Therefore, hydrogen is 
fully-ionized at the magnetosphere boundary before the 
star reaches that region]. 

Figure 1 shows the j-8 parameter space for which 
RXJ1856 can spin-down to P > 10 4 sec within 5 x 10 5 
years. Our rough interface model (7 = 5 = l,n = 3/7) 
can achieve the required spin-down for B <~ 5 x 10 15 
G (/i 33 = 5) and n m = 1 cm" 3 and B ~ 5 x 10 14 G 
(^33 = 0-5) and n m = 10 5 cm -3 (estimated density of the 
nearby molecular cloud R CrA through which RXJ1856 
may have passed (Giannini et al. 1998)). 

3.1. Thermal emission from RXJ1856 

The surface temperature of RXJ1856 is roughly con- 
sistent with predictions of standard cooling curves at its 
age - 5 x 10 5 years (Tsuruta et al. 2002). RXJ1856 
may radiate because of accretion at an incoming mass rate 
M = TrR 2 m v p m = 5x ltf R 2 m l2 v m ^ 7 n m gs" 1 . However, 
the X-ray luminosity from such an accretion rate is below 
the detection limit of Chandra and XMM-Newton obser- 
vations (Rutledgc 2001). A third heat source may be the 
continuous dissipation of the large magnetic field energy 
in the stellar crust (Heyl & Kulkarni 1998). 

Given featureless X-ray spectra, several approaches us- 
ing spectral energy distribution (SED) in the optical and 
X-ray band have been undertaken to reveal the surface 
composition (Pons et al. 2002). Light element atmospheres 
(H and He) were ruled out since they predict about two 
orders of magnitudes larger optical flux compared to the 
observed values (Pons et al. 2002). On the other hand, 



a blackbody model underpredicts optical flux by a fac- 
tor of ~ 7 (Walter & Lattimer 2002). Non-magnetized 
heavy element atmosphere models (e.g. Si-ash or Iron) 
predict correct SED over the optical and X-ray band with 
a single temperature (kT°° = 40 eV) (Walter & Lattimer 

2002) . However, they have absorption features which devi- 
ate from the featureless X-ray data (Burwitz et al. 2003). 
A similar situation exists for magnetized heavy element 
atmospheres (Rajagopal et al. 1997). 

A two-component blackbody model was proposed to ac- 
count for the multi-wavelength SED and the featureless 
X-ray spectra (Pons et al. 2002). Hot and cold blackbody 
components for X-ray and optical spectra respectively are 
consistent with both SED and featureless X-ray data (Pons 
et al. 2002; Braje & Romani 2002). However, it is hard to 
obtain blackbody-like spectra when significant atmosphere 
is present on the surface. 

At sufficiently high magnetic field strength and low tem- 
perature, a NS surface becomes very dense liquid or solid 
with almost no atmosphere above it (Ruderman 1971; Lai 
& Salpeter 1996). 

That RX J 1856 is a magnetar is consistent with the two 
temperature model (anisotropy of surface temperature dis- 
tribution caused by strong magnetic field) and the model 
for a condensed iron surface because they both require 
strong magnetic field strength on the surface. The emis- 
sivity of such a condensed matter surface in 10 15 G is about 
1/2 that of a blackbody because the strongly magnetized 

very dense (p - 560Z- 3 / 5 B% 5 gcnT 3 , with B = 10 12 B 12 
G) has essentially no cmissivity for O-mode X-rays but is 
near blackbody for E-modc ones. Then the average emis- 
sivity for the radiating surface is ~ 30-50% depending on 
photon energy and magnetic field geometry (Zane et al. 

2003) . The actual NS area should then be more than twice 
that of the apparent blackbody area so that the inferred 
NS radius becomes <: \P1Rbb (Zane et al. 2003). 

4. DISCUSSION 

Will we find other isolated NSs in the propeller phase? 
There are several NSs with discrepant supernova remnant 
and "canonical" spin-down ages (t csc i). In some cases, the 
discrepancy may be due to the fact that NSs are in a pro- 
peller phase. (Alternatively, they could have been born 
with a long spin period close to present value.) However, 
some of them have spin periods shorter than < 1 sec, so it 
is difficult for them to enter a propeller phase unless the 
ambient gas density is very large (n m 3> 1 cm -3 ). 

Detection of long pulsation periods (> 10 sec) from iso- 
lated NSs may be another indication of NSs in a propeller 
phase. Potential candidates for such isolated NSs would 
have ages between 10 3 years (enough time for spinning 
down to enter the propeller phase) and 10 6 years (still de- 
tectable thermal emission). In a magnetar, magnetic field 
decay processes can keep the magnetar more X-ray lumi- 
nous than would be the case for canonical pulsars (Heyl & 
Kulkarni 1998). Then older magnetars might be still ob- 
servable by their X-ray emission after canonical NSs have 
faded. 

Detection of spin-period variation in X-ray flux over sev- 
eral X-ray observations would explore periods > 10 sec 
for the case of RX J1856. In most cases in which RX J1856 
has spun-down to P > 10 4 sec (§3), it should at present 
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have a period of at least 10 6 sec. The chance of seeing a 
different surface aspect of the star may be small and grav- 
itational light-bending effects will smear out some spectral 
modulation. 

If an isolated NS in a propeller phase is confirmed it 
will select among propeller effect models. If a NS is in 
the propeller phase, one should be cautious in deriving NS 
ages and magnetic field strengths from P and P measure- 
ments since the conventional dipole radiation wind formu- 



lae of equation (1) and (2) are no longer valid. In such 
cases, independent determination of ages and magnetic 
field strengths (e.g. supernova remnant ages and spec- 
troscopic measurements of B-field) are needed. 

We thank J.J. Drake, R.V.E. Lovelace, J. Pringle, J.E. 
Triimper, M.H. van Kerkwijk and S. Zane for useful dis- 
cussions. 



REFERENCES 



Braje, T. M. & Romani, R. W. 2002, ApJ, 580, 1043 

Burwitz, V., Haberl, F., Neuhauser, R., Predehl, P., Triimper, J., & 

Zavlin, V. E. 2003, A&A, 399, 1109 
Davidson, K. & Ostriker, J. P. 1973, ApJ, 179, 585 
Davies, R. E., Fabian, A. C, & Pringle, J. E. 1979, MNRAS, 186, 

779 

Drake, J. J. et al. 2002, ApJ, 572, 996 

Giannini, T. et al. 1998, in ASP Conf. Ser. 132: Star Formation with 

the Infrared Space Observatory, 350 
Haberl, F., Motch, C, Buckley, D. A. H., Zickgraf, F.-J., & Pietsch, 

W. 1997, A&A, 326, 662 
Heyl, J. S. & Kulkarni, S. R. 1998, ApJ, 506, L61 
Illarionov, A. F. & Sunyaev, R. A. 1975, A&A, 39, 185 
Kaplan, D. L., van Kerkwijk, M. H., & Anderson, J. 2002, ApJ, 571, 

447 

Kouvcliotou, C. et al. 1998, Nature, 393, 235 

Lai, D. & Salpetcr, E. E. 1996, Phys. Rev. A, 53, 152 



Menou, K., Esin, A. A., Narayan, R., Garcia, M. R., Lasota, J., & 

McClintock, J. E. 1999, ApJ, 520, 276 
Pons, J. A., Walter, F. M., Lattimer, J. M., Prakash, M., Neuhauser, 

R., & An, P. 2002, ApJ, 564, 981 
Rajagopal, M., Romani, R. W., & Miller, M. C. 1997, ApJ, 479, 347 
Ransom, S. M., Gaensler, B. M., & Slane, P. O. 2002, ApJ, 570, L75 
Romanova, M. M., Toropina, O. D., Toropin, Y., & Lovelace, R. 

V. E. 2002, submitted to ApJ (astro-ph/0209548) 
Ruderman, M. A. 1971, Phys. Rev. Lett., 27, 1306 
Rutledge, R. E. 2001, ApJ, 553, 796 

Tsuruta, S., Teter, M. A., Takatsuka, T., Tatsumi, T., & Tamagaki, 

R. 2002, ApJ, 571, L143 
van Kerkwijk, M. H. & Kulkarni, S. R. 2001, A&A, 380, 221 
Walter, F. M. & Lattimer, J. 2002, ApJ, 576, L145 
Walter, F. M., Wolk, S. J., & Neuhauser, R. 1996, Nature, 379, 233 
Zane, S., Turolla, R., & Drake, J. J. 2003, astro-ph/0302197 
Zavlin, V. E., Pavlov, G. G., Sanwal, D., & Triimper, J. 2000, ApJ, 

540, L25 



Table 1 

Propeller effect model parameters. 
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A rough interface model (RIM) 
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Romanova et al. (2002) a 



a Thcir analytic model gives n = 3/5 which differs from the n = 3/4 
because they used the outflow density which depends on SlR m (see more 
details in section 2 of Romanova et al. (2002)). In their numerical model, 
n = 1.3 and = 0.8. 
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Fig. 1. — 7-5 diagram for ^33 = 5 and n m = 1 cm -3 (left) and /133 = 0.5 and n m = 10 5 cm -3 (right). The dotted region shows he 7-5 
parameter space for which RX J1856 can spin-dwon to P > 10 4 sec within 5 X 10 5 years. The asterisk refers to our favorite model (7 = 5 = 1). 
The ISM density in the vicinity of RXJ1856 has been estimated to be ~ 1 cm -3 from studies of the nearby H« nebula (van Kerkwijk & 
Kulkarni 2001). RX J1856 may have passed through the nearby molecular cloud R CrA with density > 10 5 cm~ 3 (Giannini et al. 1998). 



Table 2 

ns periods and transition ages for different pulsar spin-down mechanisms. 



Transition 


Period [sec] 


Transition time [yrs] a 


Pulsar wind — > Propeller 
Supersonic — > Subsonic 
Propeller — > Accretion 


37/^33 Tim Pm 

1.6 x lOV^n^ 6 ^ 3 
~ [10 6 ] b 


2.2xl0V33 4/3 ^ V3 ^- 2)/3 
[10 5 -10 7 ] b 
[10 5 -10 7 ] b 



Note. — p m = v m /c and fj, = 10 33 ^33 Gem 3 . ^ 33 = 1 corresponds to B = 10 15 G 
and R = 10 km. 

a Age since birth. 

b These do not usually have simple power-law forms and we show approximate ranges. 
We assumed 7 = 5=1. The parameter ranges we considered here are: B ~ 10 14 -10 15 
G, n m ~ 1-10 cm -3 and v m ~ 10-100 kms -1 . 



